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Abstract 

A two-phase model, where the plasma expansion is an isothermal one when laser irradiates and a fol- 
lowing adiabatic one after laser ends, has been proposed to predict the maximum energy of the proton 
beams induced in the ultra-intense laser-foil interactions. The hot-electron recirculation in the ultra-intense 
laser-solid interactions has been accounted in and described by the time-dependent hot-electron density con- 
tinuously in this model. The dilution effect of electron density as electrons recirculate and spread laterally 
has been considered. With our model, the scaling laws of maximum ion energy have been achieved and 
the dependence of the scaling coefficients on laser intensity, pulse duration and target thickness have been 
obtained. Some interesting results have been predicted: the adiabatic expansion is an important process of 
the ion acceleration and cannot be neglected; the whole acceleration time is about 10 - 20 times of laser 
pulse duration; the larger the laser intensity, the more sensitive the maximum ion energy to the change of 
focus radius, and so on. 
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I. INTRODUCTION 



Proton acceleration mechanisms in ultra-intense laser pulses interaction with thin solid targets 
attract more and more interest nowadays [iQjIsI]. Various models [2,4,5,6] have been presented 
to estimate the maximum energy of proton beams. However, the models given by Wilks et al. , 
Kaluza et al, Schreiber et al. and Fuchs et al. are all based on isothermal expansions of quasi- 
neutral plasmas [vj]. Robson et al. presented a two-phase temperature- varying model, where the 
hot-electron temperature first increases linearly on the pulse duration timescale and then decreases 
adiabatically with time . However, in the pulse duration, does the hot-electron temperature rise 
up linearly? That is still difficult to be validated. Generally ones assume that, when an ultra- 
intense laser pulse interacts with a solid target, the laser-produced fast electrons with a uniform 
temperature, ksTg, determined by the laser ponderomotive potential are instantly created in front 
of the target and propagate through the target coUisionlessly and then form a high energy plasma 
at the rear of the target. When the laser pulse still exists, the hot-electron temperature, ksTe = 
nieC^iy - 1), is assumed invariant due to a constant energy supply from the laser pulse, where 
7 = (1 -I- 7/1^/ 1.37)"-^ is the relativistic factor, / is the laser intensity in lO'^W/cm^, A is the laser 
wave length in /urn, is the electron mass and f/ is the pulse duration. The plasma expansion is 
an isothermal expansion. Therefore, a two-phase model different from Robson et al. is proposed 
in this article, where the plasma expansion is isothermal in the laser pulse duration and then the 
hot-electron temperature decreases as (?/?/)"*'^'^^^ [8]. 

The electron density distribution satisfies Boltzmann relationship: = neoexp{e(f>/kBTe) and 
stays a constant and time-independent in the previous models [|2t |4|, |5i lZO without hot-electron 
recirculation, where e is the elementary charge and (p is the electric potential. Therefore, with a 
little adjustment of some parameters: the acceleration time [4,|5] , the opening angle of electrons 
[I2] and electron density, n^o [2. 4, 5], Mora's result can be used to estimate the maximum energy 
of proton beams for thick targets, where the influence of hot-electron recirculation on the ion 
acceleration can be ignored. Although Robson et al. have presented a two-phase model which is 
consistent with experiments, the hot-electron recirculation is still ignored. However, Mackinnon 
et al. dgl] observed enhancement of proton acceleration by hot-electron recirculation in thin foils 
whose thickness is less than a critical value. In addition, Sentoku et al. [|lO|l predicted an equation 
to conclude the influence of electron recirculation and proved the hot-electron recirculation cannot 
be ignored in the laser-foil interactions, although they didn't propose a clear description of electron 
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recirculation and their physical picture is too simple and not clear. The assumption: the maximum 
hot-electron density for a thin foil is a constant and N times of the value for a critical target 
thickness in Sentoku et a/.'s model, is rough and unreasonable. Because there are n times of the 
electron recirculation they happen one after the other and the electron density can't jump to n 
times of the initial density. After that, Huang et al.\\l] presented a step model to describe the 
influence of the hot-electron recirculation on the laser-ion acceleration. In the step model, the 
hot-electron density rises step by step with isothermal plasma expansions. In fact, the electron 
density should rise continuously and then decrease to zero as the time tends to infinite. Therefore, 
the time-dependent hot-electron density and the electric field are necessary for the description of 
the hot-electron recirculation and the whole process of the ion acceleration. The dilution effect 
of the electron density as the electrons circulate and spread laterally should be considered but not 
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accounted in the previous models 

In Sec. HH a new two-phase model which contains three-dimensional effect (the thickness effect, 
the angular effect, which are discussed in detail by Huang et al. [|ll|l . and the dilution effect as the 
electrons circulate and spread laterally) and the hot-electron recirculation is proposed, where the 
plasma expansion is isothermal in the pulse duration and then adiabatic. The main processes of 
our model are two: first, combining the Mora's result in ref. [7] and the increase of the electron 
density in the pulse duration, with the assumption: the hot-electron temperature is a constant, the 
dependence hot-electron density, the electric field and the ion velocity on the time are obtained; 
second, with the assumption of an adiabatic expansion, the dependence of the temperature of hot 
electrons on time as proposed by Mora|l8|] has been used and then the maximum ion velocity is 
obtained easily. A most significant progress of our model is that: the time-dependent electric 
field and hot-electron density can be given easily by solving two nonlinear equations. As a result 
of the model, the duration of the time-dependent electric field at the ion front is approximately 
one to two times of the main laser pulse duration which is consistent with the result presented 
by d'Humieres, Lefebvre, Gremillet, and Malka in yfl using the particle-in-cell (PIC) simulation. 
The whole acceleration time is about 10-20 times of the laser pulse duration. And we also 
proofed that the adiabatic expansion is an important process for the ion acceleration and cannot be 
neglected. Our model can be used in the same application content as Robson's model: the target 
normal sheath acceleration of ions, however, from the above discussions, it is more reasonable and 
easily to use than their' s. 

In Sec. Uni with a proper laser absorption efficiency for thick targets, our two-phase model has 
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been compared with experiments and they are consistent as shown in Table HI The laser absorption 
stays constant with the target thickness for thick target. The laser absorption efficiency for the 
target of arbitrary thickness has been calculated by particle-in-cell (PIC) simulations [jsl], although 
there is no analytic result of that. With the laser absorption efficiency of 40% for the target of 3/im 
given by the result of PIC simulations, the comparison between our model and the experimental 
result is shown in TableHl If the laser absorption is known, for the target of arbitrary thickness, the 
maximum energy of proton beams and time-dependent electric field and electron density can all 
be obtained using our model. 

In Sec. UyI with our two-phase model, the scaling law of maximum ion energy with respect 
to laser intensity for a series of constant pulse duration has been given and discussed as shown in 
Figures [21 The dependence of maximum ion energy on target thickness, focus radius, laser pulse 
duration have been obtained. With the scaling law, some interesting results have been obtained 
and discussed in detail. Also in this section, the limits of our model have been discussed. 



II. TIME-DEPENDENT ENERGY PROTON ACCELERATION 



For convenience, the physical parameters: the time, t, the ion position, /, the ion velocity, v, the 
electron field, E, the hot-electron density, n, and the light speed, c, are normalized as Equation (1) 
in [HID . Then the normalized parameters are : r, /, u, E, h, c as shown by Equation (1) in [Hi]. 

When an ultra- intense laser pulse interacts with a solid target, the laser-produced fast electrons 
with a uniform temperature, kgTe, determined by the laser ponderomotive potential are instantly 
created in front of the target and propagate through the target and then form a high energy plasma 
at the rear of the target. Here, it is assumed that the hot electron transport is coUisionless, which 
is true for high energetic electrons, thin foils or the atomic number of the materials of the target is 
low. Hot electrons at the rear of the target can be considered to be reflected by sheath field at the 



ion front [|7, 



1211 and come back to the front of the target, because the field there is strongest. Once 



hot electrons are created, they will bounce between the ion front before the target and the ion front 
at the rear side. Since we consider the electron motion is coUisionless, the bounce of hot electrons 
will last in the whole time of the plasma expansion. When the hot electrons propagate through the 
target, the electron beam can be assumed to be in equilibrium. 
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A. Isothermal Expansion 



The hot-electron speed used is the light speed c. Here the choice of ? = is the same as that in 



the step model given by Huang and co-workers in llllll . For simplicity in the -L/c < t < ti - L/c, 
where L is the target thickness, the laser intensity is assumed to be a constant, therefore, the hot- 
electron temperature, ksTg = nigC^iy - 1), is invariant. The plasma expansion is an isothermal 
expansion. 

The fast-electron density is a function of the parameters: the acceleration time, r, the target 
thickness, L, laser intensity, /, laser focus radius, r^, the laser absorption efficiency, 77, the incidence 
angle of the laser pulse. Bin, the half-opening angle of fast electrons, 6^. The time-dependent 
electron density is assumed: 

n,(T, L, I, ri, J], Bin, 9,) = N{t, L)n,(){L, /, r^, j], Bin, 9,), 
N{t,L) = 1,t = ti = tl, 

where ri is the time when the zeroth hot-electron recirculation ends and hot electrons go forth to 
reach the rear of the target the second time, Tl = 2L/c/ V2e, here e denotes the numerical constant 
2.71828.... HgoiL, /, ri, rj, 9i„, 9^) is the hot-electron density when hot electrons return back from 
the ion front before the target and go forth to reach the rear of the target the second time and 
N{t, L) describes the increase of the maximum electron density due to electron recirculation and 
the electron generation by the laser-plasma interactions at the front of the target. 

Using Eq. (2) in [|2[], since the total number of hot electrons that propagate through the target at 
t = tL = ILIc N, = T](L)Ei/(kBT,) for t, < and A^, = r]{L)E,td{kBTj,) for h > t^, in Eq. © 
can be estimated by: 



(7 - 1)(1 + (LVrL)tan(0e))^ 



= — — . ... ^.2 '^^ - (2) 



(y-l)(l+(L*/rL)tan(02t^' ' 
where is the laser pulse focus radius, L* = L/cos{9in) is the efficient target thickness, 9i„ is the 
incidence angle of the laser pulse and 9^, ~ 17" is half-opening angle of the superathermal electrons 



which was measured by Santos et al. II12I1 . With Eq. ^ and Eq. ([3]), the three-dimensional effect 
has been accounted in through the considering of the half-opening angle of electrons, 9e ~ 
Note that the right side of Eq. @ has a factor, tijti, which is not in the right side of Eq. ©. 
For tl < tt, Ne = ri{L)Eil{kBTt,), where Ei is the energy of laser pulse. However, for ti > tt, at 
t = tl, hot electrons are still being generated by the laser-plasma interactions at the front of the 
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FIG. 1: (Color online) The time-dependent hot-electron density, acceleration field and the speed of ions 
versus t = ojpiti V2e given by the Time-Dependent Target Normal Sheath Acceleration for 9in = 30° (in 
Fig. [11(a)), Oin = 22" (in Fig. [I](b)), Og = 17" and ij = 40%. In Figure [na),the laser pulse parameters are 
/ = 1.0 X lO^'^W/cm^ A = 790nm, rt = 2.5//m, L = 30/im and ti = 150fs. In Figure^b), the laser pulse 
parameters are 7 = 1.0 x lO^^W/cm^, A = 800nm, ri = 2.5/im, L = 3/im and ti = lOOfs. 



target, and the number of hot electrons which propagate through the target is a part of the total 
number, Ng = rjiQEiti/ikBTeti). Therefore, Equations ^ and ([3]) are obtained. When » L 
and tan{6in) «; 1, (1 + (L* / rL)tan{6i„))^ « 1, the angular effect can be neglected. Therefore, the 
influence of t]{L) and electron recirculation become dominated for thin targets. For example, for 
/ = 3 X 10^°W/cm^, A = 1.053nm and tj = 500fs [Il3n . the temperature of hot electrons is about 
5.5MeV. For //(L) = 50%, = 5pim, with Eq. ([3]), the electron density is about 4.8 x lO^^cm'^ 
and it is about the value measured by x-ray in jl^ . Since the laser pulse duration is ti, with similar 
discussion in Jill], the critical target thickness for the hot-electron recirculation is: = 0.5cf/.. 

With reference to the discussion and method given by Huang and co-workers in [11], the re- 
lationship between the ion velocity at the ion front and the electron density can be described by 
Equation (12) in lllin . With that equation, the ion velocity is decided by N{t}. Therefore, the 
solution of N{t) is a key point. Although it has been given by Huang et al. jllll with a simple 
model, it is rough for three reasons: (1) it is discrete, however the actual electron density change 
continuously; (2) the electron density deceases as electrons recirculate and spread laterally and 
this dilute effect was not counted in the Step Model ill ill : (3) the turning point of hot electrons at 
front of the target is not static but moving with the expansion of the plasma too. In this paper, a 
more actual and valuable method will be proposed to calculate a continuous solution of N{t) as 
follows. 
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Assuming that: the velocity at the ion front before the target is the same as that at the rear of the 
target and the two turning points for the electron recirculation are the ion front before the target 
and at the rear respectively, if hot electrons satisfy uniform distribution in the bulk from the ion 
front before the target to the ion front at the rear, N{t) is decided by 

XT-LI c 

^(T) = y- -y— ^Fe,3-D, T < Ti, (4) 

£r "fdT 2/(T) + L 

where / represents the generation rate of hot electrons in the interaction of laser pulses with the 
plasma at the front of the target, Fe,3-D = ^i^d = tm(6)/ cos(6i„). f depends on the 

absorption mechanisms of laser pulses and decides the density of hot electrons. The factor, F^^-z), 
corresponds to the decrease of on-axis density as hot electrons circulate and spread laterally with 
an given opening angle, 9^. And also it reflects three-dimensional effect on the ion acceleration. 
The special integrating limits are because there is an interval before the electrons generated by the 
laser pulse at the front of the target come to the rear. With the assumption: / = / for t e [0, t/], 
Eq. (HI) can be simplified to be: 

T L + 21(tl) ^ ... 

N(t) = -Fe,3-D, T < Ti, (5) 

tl 2Z(t) + L 

B. Adiabatic Expansion 

When t > ti, the laser pulse has gone and the acceleration field at the ion front decreases quickly 
for two reasons. First, the temperature of hot electrons decreases with time as shown by[18|]: 

Te {TlTir''^'l^\ (6) 



For the nonrelativistic case, 7=1, with Eg. T . 



work of adiabatic expansion into a vacuumlll4 



15 



oc t ^ which is consistent with all the previous 



16h . For the ultra-relativistic case, y — > +oo, 



with Eq. (|6l), oc which is the same with Mora's results[|81]. After the laser pulse vanishes, the 
ion front does not stop and the electron bulk still increases. Therefore, the electron density, N{t), 
decreases as given by: 

T, L + 21{tl) ^ ^ 

N{t) = -Fg^^.o, r > Ti. (7) 

tl 2Z(t) + L 

Eq. © and Eq. ^ are all nonlinear differential equations and have no analytic solutions. 
However, the numerical results can be obtained by computer with iterative method. The initial 
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N{t) is given by the solution of the Eq. ([5]) and Eq. ^ in which Fq^^-d = 1- As an example, the 
solutions of a thin foil and a thick solid target have been given by Figure[T] Figure[TJa) corresponds 
a thick target of 30/im and the hot electron recirculation can be ignored. Figure [IJb) corresponds 
a thin foil of 3/im and the maximum value of N{t) is about 1.5, which is lower than that given by 
Huang and co-workers in and Sentoku et al. in [fioll and reflects the three-dimensional effect. 



From Fig. [B some interesting results can be obtained: 

(1) The whole accelerate time of ions is about 10 - 20 times of the laser pulse duration. After 
that, the separating field is close to zero and the acceleration ends. 

(2) The electron density and the electric field reach their maximum value at the time f = f/ as 
expected by our discussion and the gain energy of ions in the process of the isothermal expansion is 
approximately a quarter of the finally energy. Therefore, the adiabatic expansion is also important 
for the ion acceleration although the electron density and electric field decrease in this process. 

(3) The influence of the hot-electron recirculation on the ion acceleration for thin foils is obvi- 
ous. The maximum ion energy for thin foils is larger than that for thick targets. 

With solutions of N{t), the time-dependent electric field and the ion velocity at the ion front 
can be obtained. Therefore, for the target of arbitrary thickness, the maximum energy of proton 
beams can be achieved if the laser absorption efficiency is known. 

ni. COMPARISON WITH EXPERIMENTS 

Our time-dependent model is compared with experiments, the results are listed in Table HI 
For example, for / = 1 x lO^^W/cm^, A = 0.8nm and ti = lOOfs, the critical target thickness is 
about \5fim, according to Machinnon et al. ^ , E,„axiL = 30//m) = 6.2MeV. With the simulation 
results (Figure 12 in [3]), the laser absorption stays constant of about 35% 50% with the 
target thickness for thick target, L ^ l/zm. The laser absorption changes with target thickness, the 
contrast ratio between the main pulse and prepulse and the prepulse duration ||2Llll|] for L ^ 1/im. 
For diff'erent target thickness, the permeation of laser pulse is difl^erent. For difi'erent contrast ratio 
and prepulse duration, the scale length of the preplasma is different, which induces different laser 
absorption mechanism. Therefore, the laser absorption efficiency, r](L), is different and difficult to 
be assured. Different rj{L) corresponds to different electron density, UeQ. The plasma frequency and 
acceleration parameters depend on hm. After all, the changing law of ri(L) with L for L < l/um is 
quite important for the proton acceleration and still a challenge. Without ri{L), our model can not 
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TABLE I: This is a comparison between our two-phase model and some experiments for 0i„ = 30" in ref. 
m, Oin = 22" in ref. [9] and 9^ = \T\12]. 



I 


A 


ti 


rt 


L 






Emax from experiments 


Emax from our Model 


(lO^^W/cm^) 


(/im) 


(fs) 


(yum) 


(yum) 


(%) 


(lO^O/cm^) 


or PIC (MeV) 


(MeV) 


10 


0.79 


150 


2.5 


30 


40 


0.33 


1.2 + 0.3 [2] 


1.1 


10 


0.79 


150 


2.5 


20 


40 


0.81 


2.0±0.3[2] 


2.0 


13 


0.79 


150 


2.5 


30 


40 


0.36 


1.5 ±0.3[2] 


1.4 


15 


0.79 


150 


2.5 


30 


40 


0.37 


1.7 ±0.3[2] 


1.6 


100 


0.8 


100 


2.5 


3 


40 


14.2 


22 - 24[9] 


22.6 


100 


0.8 


100 


2.5 


6 


40 


8.6 


17 - 19[9J 


17.3 


100 


0.8 


100 


2.5 


10 


40 


5.13 


11 - iim 


13.2 


100 


0.8 


100 


2.5 


25 


40 


0.897 


6 - 7[9] 


5.0 



been compared with experiments for L < lyum. However, the laser absorption efficiency for the 
target of arbitrary thickness has been calculated by particle-in-cell (PIC) simulations [sl], although 
there is no analytic result of that. ^ 

With the simulation results in ref. [3], the small target thickness will lead to reduced absorption 
if the target deconfines rapidly and becomes transparent before the end of the laser pulse - but the 
characteristic velocity for this is the sound speed, not the speed of light. If the critical thickness 
for recirculation is L^, the critical thickness for modified absorption should be much smaller than 
Lc. Therefore, for thin foils of the thickness, L > l//m, the laser absorption efficiency keeps a 
constant about 35 - 50% approximately with the thickness, L. For the target of 3yum, rj ^ 40% and 
din = 17°, the maximum proton energy is 22.6MeV estimated by our model which is consistent 
with the experimental data, 22 - 24MeV. The time-dependent electron density, the acceleration 
field and the ion speed are shown by Figure [IJb). The hot-electron density increases from and 
reaches the maximum value 1.5 at the time Therefore, for L = 3//m and = l5/um, the hot- 
electron recirculation does exist and N(t) is up to about 1.5 but not 5 as shown by Sentoku et al. 

o 

HIOH . After ti, the electron density decreases quickly to half at about 2.2?/. The duration of the 
hot-electron density and field are 2 - 5ti. 

The maximum value of N{t) is smaller than or equal to 1 as shown by Figure [IJ Figure [T] shows 
that: the maximum N{t) is about 0.79 which is smaller than 1, therefore there is no hot-electron 
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recirculation phenomena for L> L^ ; as the time goes to infinite, the velocity of protons is finite and 
the maximum energy is about 1.1 MeV while the experimental data is 1.2 + 0.3MeV; the duration 
of the hot-electron density is about 3.6?;. Therefore, the duration of the field at the ion front is 
about 2ti, which is consistent with the simulation resultyfl. 



IV. SCALING LAW AND DISCUSSION 



The laser intensity in our model is assumed to be a constant value in the pulse duration. Under 
this assumption and for a fixed laser energy, the dependence of maximum ion energy on the laser 
pulse duration is easy to be obtained. There is an optimum pulse duration for the target normal 
sheath acceleration of ions if the laser energy, focus radius and absorption efficiency sustain con- 
stants. It is a conflict of large acceleration gradient and long efficient acceleration time. For long 
pulse duration, the intensity will be low and the acceleration field will be low. For a high intensity, 
the efficient acceleration time will be short. Therefore, there is an optimum pulse duration in the 
ion acceleration. 

For diff'erent focus radius and target thickness, the dependence of maximum ion energy, E^axh 
on pulse duration can be also obtained easily with our model. The results may be diff"erent since 
the plasma density changes with and L as shown by Eq. Q and Eq. ([3]). The wave length will 
not influence the dependence of E,„ax,i on pulse duration. 

For a series of given pulse duration, the dependence of E,nax,i on laser intensity has also been 
obtained and the scaling law is given by: 

{eXp(ai)(/iol8w/cm2^^m) ''Ao'S W/cm'^'^pm ^ 6.4, 
eXp(a2)(/l0I8w/cm2^pm)''^ Aol8W/cm2^pm ^ ^.4, 

where ai, a2, Z?2 are all coefficients and shown by Fig. [2l 

With Eq. dHl) and Fig. |2l two important results can be obtained. First, the scaling law is 
different from the previous results /^^^. The index bi and b2 depend on the laser pulse duration and 
decrease with pulse duration. It shows the adiabatic expansion of plasmas is also a very important 
acceleration process and should not be neglected, although they[|2|,|5D can consist with experiments 
considering the isothermal expansion only through adjusting of the parameters: 77, the acceleration 
time tacc, the plasma density n^, the opening angle of electrons and so on. Second, the influence of 
hot-electron recirculation on the ion acceleration can be shown by Fig. [21 approximately. 
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FIG. 2: (Color online) The coefficients in the scaling law given by Eq. [8] versus laser intensity, /. Laser 
absorption efficiency is assumed 40%, since the target thickness is large enough. Here, focus radius is 
2.5/zm. 

With our model, the dependence of maximum ion energy on target thickness is given by: 



L-Lo 

Emax,i =yo+A exp( — ), 



(9) 



for T] = 40%, = l.S/um, = 22", ti = lOOfs and 6, = 17°, where yo. A, Lq, L, are all 
coefficients. The dependence of them on laser intensity have been shown in Fig. [21 Eq. dH) shows 
the maximum ion energy decreases with the target thickness in the negative exponential form for 
a fixed laser absorption and a given laser intensity. However, the laser absorption tends to zero 
as L ^ 0, therefore, the maximum ion energy tends to zero in fact. From Fig. [31 maximum 
ion energy increases with the laser intensity since the coefficients yQ,A all increase with the laser 
intensity. 

In fact, focus radius influences the electron density and then the ion acceleration. For L = 
10/im, i] = 40%, 9in = 22", 9e = 17", using our model, the effect of focus radius satisfies: 



Emax,i — ■A2 + 



(10) 



where Ai, A2, and p are all coefficients and change with laser intensity and shown by Fig. [H 
With Eq. ([TOI ). some interesting results can be achieved: 

(1) The smaller laser intensity, the larger the critical value and the index p. Therefore, for 
ri » ro and (r^/ro)'' » 1, £,„a.v,,- « -(rL/ro)'''. The derivation of E,„aj,/. dE,riax,i/dirL/ro) = 
pirt/roy^^, reflects the rate of change of E,„ax,i to the focus radius. The rate of change is positive 
but decreases with focus radius. For larger laser intensity, the larger the rate of change. 
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FIG. 3: (Color online) The coefficients of the scaling law of maximum ion energy with respect to target 
thickness given by Eq. ^ versus laser intensity for = 2.5yum, rj = 40%, ti = lOOfs and A = 0.8/im. 




100 200 300 400 100 200 300 400 

18 2 

laser intensity (10 W/cnn ) 



FIG. 4: (Color online) The coefficients of the scaling law of maximum ion energy with respect to focus 
radius given by Eq. (flOl ) versus laser intensity for L = lOjum, rj = 40%, ti = lOOfs and A = 0.8/im. 



(2) Oppositely, for rt ro and {ri/roy «; 1, Emaxj * constant, which shows the influence of 
ri on the maximum ion energy can be ignored in this case. For larger laser intensity, the smaller 
critical focus radius. Therefore, the larger laser intensity, the more sensitive the ion acceleration 
to the change of focus radius. 

The influence of opening angle of hot electrons on maximum ion energy has been discussed 
in detail by Huang and co-workers in ref. [|ll|l . However, the dilution of the electron density as 
electrons circulate and spread laterally was not contained there. This effect is considered here with 
the factor, Fg^^-o in Equations ©-(IT]). Therefore, the maximum value of electron density here is 
1.5 (in Fig. \V[b)) while it is about 5 in ref. [H] and ref. Iioll . For the target of arbitrary thickness, 
the maximum energy of ions heated by target normal sheath acceleration (TNSA) can be obtained 
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by this model if the absorption efficiency of laser pulse is given. 

Here we will discuss the limits of our model. The prepulse is not considered in our model and 
the contrast is assumed large about 10^ which can be achieved in lots of experiments. However, 
the exist of a prepulse would generate a preplasma and the preplasma size, the scaling length of 
the preplasma, will most influence the mechanisms of laser absorption and then the temperature 
of hot electrons. Different laser absorption mechanism results in different generation rate of hot 
electrons, /, and different hot-electron temperature, Te. No matter what the mechanisms are, the 
generation of hot electrons is cumulative and the assumption: f = f causes little error relative 
to that caused by the measurement in experiments. Whatever the temperature of hot electrons is, 
our model is still in use with the actual temperature instead of the value, mc^{y - 1). The laser 
intensity in our model is assumed to be a constant value in the pulse duration. In fact, the intensity 
is changing with time and the distribution is about Gaussian distribution. However, the actual 
intensity distribution with respect to time and position when laser pulse is acting on a target is 
quite difficult to be measured in real time. Since we do not consider the time distribution of laser 
intensity, we can not give an estimation of the error. In the next paper, we will consider that case. 

V. CONCLUSION 

In conclusion, the time-dependent isothermal expansion and adiabatic expansion for the target 
normal sheath proton acceleration is discussed. A two-phase model and a new scaling law of the 
maximum energy of proton beams have been proposed. The influence of the hot-electron recir- 
culation on the ion acceleration has been accounted in. For L > L^, the hot-electron recirculation 
can be ignored. But for L < Lc, the hot-electron recirculation exists and enhances as the target 
thickness decreases. The results given by our model have been compared with experiments and 
shown in TableHl The dependence of maximum ion energy on target thickness, focus radius, laser 
pulse duration have been obtained and shown by equations ©, dH), (flOl) . and Figures El S |4] and 
so on. At last, the application and limits of our model has been discussed. 

However, for thin foils, the laser absorption efficiency is an important parameter for our model 
and is still a challenge for this problem. The generation rate of hot electrons in the interaction of 
laser pulses with the plasma at the front of the target, /, is also a challenge. An interesting work 
that may be considered nest is the time-dependent laser pulse intensity in order to optimize our 
model further more. 
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